Search Results

High-resolution planar laser-induced fluorescence (PLIF) measurements were performed on the scalar field in an optical engine. The measurements were of sufficient resolution to fully resolve all of the length scales of the flow field through the full cycle. The scalar dissipation spectrum was calculated, and by fitting the results to a model turbulent spectrum the Batchelor scale of the turbulent flow was estimated. The scalar inhomogeneity was introduced by a low-momentum gas jet injection. A consistent trend was observed in all data; the Batchelor scale showed a minimum value at top dead center (TDC) and was nearly symmetric about TDC. Increasing the engine speed resulted in a decrease of the Batchelor scale, and the presence of a shroud on the intake valve, which increased the turbulence intensity, also reduced the Batchelor scale. The effect of the shrouded valve was less significant compared to the effect of engine speed.

The paper presents the development of a novel approach to the solution of detailed chemistry in internal combustion engine simulations, which relies on the analytical computation of the ordinary differential equations (ODE) system Jacobian matrix in sparse form. Arbitrary reaction behaviors in either Arrhenius, third-body or fall-off formulations can be considered, and thermodynamic gas-phase mixture properties are evaluated according to the well-established 7-coefficient JANAF polynomial form. The current work presents a full validation of the new chemistry solver when coupled to the KIVA-4 code, through modeling of a single cylinder Caterpillar 3401 heavy-duty engine, running in two-stage combustion mode.

An investigation of high speed direct injection (DI) compression ignition (CI) engine combustion fueled with gasoline injected using a triple-pulse strategy in the low temperature combustion (LTC) regime is presented. This work aims to extend the operation ranges for a light-duty diesel engine, operating on gasoline, that have been identified in previous work via extended controllability of the injection process. The single-cylinder engine (SCE) was operated at full load (16 bar IMEP, 2500 rev/min) and computational simulations of the in-cylinder processes were performed using a multi-dimensional CFD code, KIVA-ERC-Chemkin, that features improved sub-models and the Chemkin library. The oxidation chemistry of the fuel was calculated using a reduced mechanism for primary reference fuel combustion chosen to match ignition characteristics of the gasoline fuel used for the SCE experiments.

Many recent studies have shown that the Reactivity Controlled Compression Ignition (RCCI) combustion strategy can achieve high efficiency with low emissions. However, it has also been revealed that RCCI combustion is difficult at high loads due to its premixed nature. To operate at moderate to high loads with gasoline/diesel dual fuel, high amounts of EGR or an ultra low compression ratio have shown to be required. Considering that both of these approaches inherently lower thermodynamic efficiency, in this study natural gas was utilized as a replacement for gasoline as the low-reactivity fuel. Due to the lower reactivity (i.e., higher octane number) of natural gas compared to gasoline, it was hypothesized to be a better fuel for RCCI combustion, in which a large reactivity gradient between the two fuels is beneficial in controlling the maximum pressure rise rate.

The focus of this paper is to discuss the modeling and control of intake plenum pressure on the Powertrain Control Research Laboratory's (PCRL) Single-Cylinder Engine (SCE) transient test system using a patented device known as the Intake Air Simulator (IAS), which dynamically controls the intake plenum pressure, and, subsequently, the instantaneous airflow into the cylinder. The IAS exists as just one of many devices that the PCRL uses to control the dynamic boundary conditions of its SCE transient test system to make it “think” and operate as though it were part of a Multi-Cylinder Engine (MCE) test system. The model described in this paper will be used to design a second generation of this device that utilizes both continuously and discretely actuating valves working in parallel.

The simulation of combustion chemistry in internal combustion engines is challenging due to the need to include detailed reaction mechanisms to describe the engine physics. Computational times needed for coupling full chemistry to CFD simulations are still too computationally demanding, even when distributed computer systems are exploited. For these reasons the present paper proposes a time scale separation approach for the integration of the chemistry differential equations and applies it in an engine CFD code. The time scale separation is achieved through the estimation of a characteristic time for each of the species and the introduction of a sampling timestep, wherein the chemistry is subcycled during the overall integration. This allows explicit integration of the system to be carried out, and the step size is governed by tolerance requirements.

The paper describes a methodology to co-simulate, with high fidelity, simultaneously and in one computational framework, all of the main vehicle subsystems for improved engineering design. The co-simulation based approach integrates in MATLAB/Simulink a physics-based tire model with high fidelity vehicle dynamics model and an accurate powertrain model allowing insights into 1) how the dynamics of a vehicle affect fuel consumption, quality of emission and vehicle control strategies and 2) how the choice of powertrain systems influence the dynamics of the vehicle; for instance how the variations in drive shaft torque affects vehicle handling, the maximum achievable acceleration of the vehicle, etc. The goal of developing this co-simulation framework is to capture the interaction between powertrain and rest of the vehicle in order to better predict, through simulation, the overall dynamics of the vehicle.

Two different types of mesh used for diesel combustion with the KIVA-4 code are compared. One is a well established conventional KIVA-3 type polar mesh. The other is a non-polar mesh with uniform size throughout the piston bowl so as to reduce the number of cells and to improve the quality of the cell shapes around the cylinder axis which can contain many fuel droplets that affect prediction accuracy and the computational time. This mesh is specialized for the KIVA-4 code which employs an unstructured mesh. To prevent dramatic changes in spray penetration caused by the difference in cell size between the two types of mesh, a recently developed spray model which reduces mesh dependency of the droplet behavior has been implemented. For the ignition and combustion models, the Shell model and characteristic time combustion (CTC) model are employed.

The ability to make fully resolved turbulent scalar field measurements has been demonstrated in an internal combustion engine using one-dimensional fluorobenzene fluorescence measurements. Data were acquired during the intake stroke in a motored engine that had been modified such that each intake valve was fed independently, and one of the two intake streams was seeded with the fluorescent tracer. The scalar energy spectra displayed a significant inertial subrange that had a −5/3 wavenumber power dependence. The scalar dissipation spectra were found to extend in the high-wavenumber regime, to where the magnitude was more than two decades below the peak value, which indicates that for all practical purposes the measurements faithfully represent all of the scalar dissipation in the flow.

A study of partially premixed combustion (PPC) with non-oxygenated 91 pump octane number1 (PON) commercially available gasoline was performed using a heavy-duty (HD) compression-ignition (CI) 2.44 l Caterpillar 3401E single-cylinder oil test engine (SCOTE). The experimental conditions selected were a net indicated mean effective pressure (IMEP) of 11.5 bar, an engine speed of 1300 rev/min, an intake temperature of 40°C with intake and exhaust pressures of 200 and 207 kPa, respectively. The baseline case for all studies presented had 0% exhaust gas recirculation (EGR), used a dual injection strategy a -137 deg ATDC pilot SOI and a -6 deg ATDC main start-of-injection (SOI) timing with a 30/70% pilot/main fuel split for a total of 5.3 kg/h fueling (equating to approximately 50% load). Combustion and emissions characteristics were explored relative to the baseline case by sweeping main and pilot SOI timings, injection split fuel percentage, intake pressure, load and EGR levels.

A continuous multi-component fuel evaporation model has been integrated with an improved G-equation combustion and detailed chemical kinetics model. The integrated code has been successfully used to simulate a gasoline direct injection engine. In the multi-component fuel model, the theory of continuous thermodynamics is used to model the properties and composition of multi-component fuels such as gasoline. In the improved G-equation combustion model a flamelet approach based on the G-equation is used that considers multi-component fuel effects. To precisely calculate the local and instantaneous residual which has a great effect on the laminar flame speed, a “transport equation residual” model is used. A Damkohler number criterion is used to determine the combustion mode in flame containing cells.

A transport equation residual model incorporating refined G-equation and detailed chemical kinetics combustion models has been developed and implemented in the ERC KIVA-3V release2 code for Gasoline Direct Injection (GDI) engine simulations for better predictions of flame propagation. In the transport equation residual model a fictitious species concept is introduced to account for the residual gases in the cylinder, which have a great effect on the laminar flame speed. The residual gases include CO2, H2O and N2 remaining from the previous engine cycle or introduced using EGR. This pseudo species is described by a transport equation. The transport equation residual model differentiates between CO2 and H2O from the previous engine cycle or EGR and that which is from the combustion products of the current engine cycle.

The effect of the ring pack storage mechanism on the hydrocarbon (HC) emissions from an air-cooled utility engine has been studied using a simplified ring pack model. Tests were performed for a range of engine load, two engine speeds, varied air-fuel ratio and with a fixed ignition timing using a homogeneous, pre-vaporized fuel mixture system. The integrated mass of HC leaving the crevices from the end of combustion (the crank angle that the cumulative burn fraction reached 90%) to exhaust valve closing was taken to represent the potential contribution of the ring pack to the overall HC emissions; post-oxidation in the cylinder will consume some of this mass. Time-resolved exhaust HC concentration measurements were also performed, and the instantaneous exhaust HC mass flow rate was determined using the measured exhaust and cylinder pressure.

This paper investigates flow and combustion in a full-cycle simulation of a four-stroke, three-valve HCCI engine by visualizing the flow with pathlines. Pathlines trace massless particles in a transient flow field. In addition to visualization, pathlines are used here to trace the history, or evolution, of flow fields and species. In this study evolution is followed from the intake port through combustion. Pathline analysis follows packets of intake charge in time and space from induction through combustion. The local scalar fields traversed by the individual packets in terms of velocity magnitude, turbulence, species concentration and temperatures are extracted from the simulation results. The results show how the intake event establishes local chemical and thermal environments in-cylinder and how the species respond (chemically react) to the local field.

The lift-off length plays a significant role in spray combustion as it influences the air entrainment upstream of the lift-off location and hence the soot formation. Accurate prediction of lift-off length thus becomes a prerequisite for accurate soot prediction in lifted flames. In the present study, KIVA-3v coupled with CHEMKIN, as developed at the Engine Research Center (ERC), is used as the CFD model. Experimental data from the Sandia National Labs. is used for validating the model predictions of n-heptane lift-off lengths and soot formation details in a constant volume combustion chamber. It is seen that the model predictions, in terms of lift-off length and soot mass, agree well with the experimental results for low ambient density (14.8 kg/m3) cases with different EGR rates (21% O2 - 8% O2). However, for high density cases (30 kg/m3) with different EGR rates (15% O2 - 8% O2) disagreements were found.

Multi-zone CFD simulations with detailed kinetics were used to model iso-octane HCCI experiments performed on a single-cylinder research engine. The modeling goals were to validate the method (multi-zone combustion modeling) and the reaction mechanism (LLNL 857 species iso-octane) by comparing model results to detailed exhaust speciation data, which was obtained with gas chromatography. The model is compared to experiments run at 1200 RPM and 1.35 bar boost pressure over an equivalence ratio range from 0.08 to 0.28. Fuel was introduced far upstream to ensure fuel and air homogeneity prior to entering the 13.8:1 compression ratio, shallow-bowl combustion chamber of this 4-stroke engine. The CFD grid incorporated a very detailed representation of the crevices, including the top-land ring crevice and head-gasket crevice. The ring crevice is resolved all the way into the ring pocket volume. The detailed grid was required to capture regions where emission species are formed and retained.

Transient engine operations are modeled and simulated with a 1-D code (GT Power) using heat release and emission data computed by a 3-D CFD code (Kiva3). During each iteration step of a transient engine simulation, the 1-D code utilizes the 3-D data to interpolate the values for heat release and emissions. The 3-D CFD computations were performed for the compression and combustion stroke of strategically chosen engine operating points considering engine speed, torque and excess air. The 3-D inlet conditions were obtained from the 1-D code, which utilized 3-D heat release data from the previous 1-D unsteady computations. In most cases, only two different sets of 3-D input data are needed to interpolate the transient phase between two engine operating points. This keeps the computation time at a reasonable level. The results are demonstrated on the load response of a generator which is driven by a medium-speed diesel engine.

Optimal fuel injection strategies are obtained with a micro-genetic algorithm and an adaptive gradient method for a nonroad, medium-speed DI diesel engine equipped with a multi-orifice, asynchronous fuel injection system. The gradient optimization utilizes a fast-converging backtracking algorithm and an adaptive cost function which is based on the penalty method, where the penalty coefficient is increased after every line search. The micro-genetic algorithm uses parameter combinations of the best two individuals in each generation until a local convergence is achieved, and then generates a random population to continue the global search. The optimizations have been performed for a two pulse fuel injection strategy where the optimization parameters are the injection timings and the nozzle orifice diameters.

Fundamental simulations in a quiescent cell under adiabatic conditions were made to understand the effect of temperature, equivalence ratio and the components of the recirculated exhaust gas, viz., CO2 and H2O, on the combustion of n-Heptane. Simulations were made in single phase in which evaporated n-Heptane was uniformly distributed in the domain. Computations were made for two different temperatures and four different EGR levels. CO2 or H2O or N2was used as EGR. It was found that the initiation of the main combustion process was primarily determined by two competing factors, i.e., the amount of initial OH concentration in the domain and the specific heat of the mixture. Further, initial OH concentration can be controlled by the manipulating the ambient temperature in the domain, and the specific heat capacity of the mixture via the mixture composition. In addition to these, the pre combustion and the subsequent post combustion can also be controlled via the equivalence ratio.

The use of virtual prototyping early in the design stage of a product has gained popularity due to reduced cost and time to market. The state of the art in vehicle simulation has reached a level where full vehicles are analyzed through simulation but major difficulties continue to be present in interfacing the vehicle model with accurate powertrain models and in developing adequate formulations for the contact between tire and terrain (specifically, scenarios such as tire sliding on ice and rolling on sand or other very deformable surfaces). The proposed work focuses on developing a ground vehicle simulation capability by combining several third party packages for vehicle simulation, tire simulation, and powertrain simulation. The long-term goal of this project consists in promoting the Digital Car idea through the development of a reliable and robust simulation capability that will enhance the understanding and control of off-road vehicle performance.